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Abstract 



X. 

^ . Using data collected with the L3 detector at LEP from 1992 to 1995 on the Z 

Qh! peak, we determine the branching fractions of the r lepton into one, three and five 

(-H I charged particles to be: 
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The first uncertainties are statistical and the second systematic. The accuracy of 
these measurements alone is similar to that of the current world average. 
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= 85.274 ±0.105 ±0.073%, 
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= 14.556 ±0.105 ±0.076%, 
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= 0.170 ±0.022 ±0.026%. 
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Introduction 

Measurements of the topological branching fractions of the r lepton and the sum of measure- 
ments of the exclusive branching fractions were previously inconsistent. Solving this "one-prong 
puzzle" motivated many precise determinations of the exclusive r branching fractions at the 
permille level jT] but only a few less precise determinations of the topological branching frac- 
tions |2] have been performed. 

In this letter, we present a new measurement of the topological branching fractions using 
data collected by the L3 detector at LEP on the Z resonance. These results supersede those of 
our previous publication j3]. Here we follow the convention that tracks stemming from neutral 
kaon decays are not accounted for in the topology. Another measurement of the topological r 
branching fractions was recently reported in Reference H| 

The measurement entails a selection of e^e~ -^ T'^T~{-f) events followed by an event topol- 
ogy reconstruction, which must be precisely understood. The reconstructed topology is in- 
fluenced by photon conversions, subdetector inefficiencies and resolution limitations. Detailed 
studies of these effects are performed in order to control the systematic uncertainties to the 
level of the statistical uncertainties. 

Data and Monte Carlo samples 

The data used were collected with the L3 detector [5^ at LEP from 1992 to 1995 on the Z peak, 
corresponding to an integrated luminosity of 92.6 pb~^. The most crucial subdetectors for this 
analysis are: the central tracking system consisting of a silicon microvertex detector (SMD), 
a time expansion chamber (TEC) and proportional chambers measuring the Z coordinate, 
the electromagnetic calorimeter composed of Bismuth Germanium Oxide (BGO) crystals, the 
hadron calorimeter (HCAL) and the muon spectrometer. Detailed studies of the efficiencies of 
these subdetectors using control samples are performed, yielding precise determinations of the 
efficiencies. 

For efficiency studies, e"'"e'^ -^ r+r^(7) events are generated with the KORALZ Monte 
Carlo generator j^. Background estimations are performed using the following Monte Carlo 
generators: KORALZ for e+e" -^ ^+/U"(7); BHAGENE for e+e" -^ e+e-{-f); JETSET 
for e"'"e~ — > qq(7); DIAG36 P for e"'"e~ — > e'^e~i~^i~, where i = e, n, or r. The Monte 
Carlo events are simulated in the L3 detector using the GEANT program ^Oj, which takes 
into account the effects of energy loss, multiple scattering and showering. Furthermore, time 
dependent detector inefficiencies are considered. These events are reconstructed with the same 
program as the one used for the data. 

Subdetector efficiencies and calibrations 

Efficiency studies of the subdetectors are done separately for each year of data taking. As the 
year-by-year efficiency variations are small, average values are given in the following. 

The efficiency of the TEC to measure a track is studied using data samples of Bhabha and 
dimuon events and muons originating from r decays. The Bhabha and dimuon samples are 
selected by requiring two energy deposits in the BGO or two tracks in the muon spectrometer 
of about the beam energy and back-to-back topology. The muons in r decays are identified as 
tracks in the muon chambers pointing to the interaction region. In addition, the energy deposits 



in the BGO and HCAL must be consistent with the expectation for a minimum ionising particle 
(MIP). 

A track in the TEC must have at least 25 out of the 62 possible hits, one or more hits in the 
innermost part of the chamber and to span over more than 40 anode wires radially. Its transverse 
momentum, p^, must be larger than 2 GeV. After rejecting tracks in the low resolution region 
adjacent to the anode, the track finding efficiency is found to be about 96%, almost independent 
of the track momentum. The double track resolution of the TEC is determined from data ^T] 
to be about 500/im and is modeled correspondingly in the detector simulation. As a cross check, 
the distributions of the azimuthal angle between two adjacent tracks A0, from data and Monte 
Carlo are compared in Figure ^ for small values of A(f) and high momentum tracks from 3-prong 
r decays. For A0 larger than 0.005 rad excellent agreement is found. The small discrepancy 
below 0.005 rad is taken into account as a systematic uncertainty. 

The efficiency of the BGO to detect an electromagnetically showering particle is determined 
using Bhabha events to be about 99.5%. This efficiency is found to be almost independent 
of the shower energy from studies using e"'"e~ —>■ e"'"e~e"'"e~ events. In order to estimate the 
efficiency of the BGO to detect a MIP, r decays into muons are used. A track in the muon 
spectrometer is required, which points to the interaction region and matches an energy deposit 
in the HCAL that corresponds to a MIP. From these muons 97% induce a signal in the BGO. 

The same technique was used to estimate the HCAL and muon spectrometer efficiencies. 
Using muons with a track in the TEC, a MIP signal in the BGO and a matched muon spec- 
trometer track, the efficiency of the HCAL to detect such a particle is about 89%. The muon 
spectrometer efficiency is found to be 74% using r decays with a track in the TEC and a MIP 
signature in the BGO and the HCAL. 

The subdetector efficiencies obtained from each year are used to correct the Monte Carlo 
simulation of the detector response for the e+e" — > r+r~(7) and background processes. The 
energy scales of the subdetectors are calibrated using control data samples ^21- The momentum 
scale of the central tracker is verified to 0.5% from 1 to 45 GeV. The BGO and the muon 
spectrometer scale uncertainties are 0.5% at low energies and 0.05% at high energy. The scale 
uncertainty of the HCAL is 1%. 

Study of photon conversions 

Photon conversions occurring in the material inside the TEC may cause additional tracks 
and are studied on data and Monte Carlo for each year independently. A loose selection of 
e"'"e~ — > r+r^(7) events is made, requiring two low multiplicity jets and the cosine of the polar 
angle of the event thrust axis | cosOthmstl < 0.7. Radiative photons or photons from tt" decays 
can convert in the detector. The tracks from the conversion point either to their corresponding 
cluster in the BGO calorimeter, or to a coalescent cluster including the energy of the two 
conversion tracks or of the 2 photons in the case of a 7r° decay. Therefore, the pt measured in 
the central tracker must be smaller than the transverse energy observed in the electromagnetic 
calorimeter. 

In the case that only one track reaches the TEC, its transverse momentum must be less than 
4 GeV. When both tracks are reconstructed, the square of their invariant mass must be less 
than 0.005 GeV^. Taking track pairs which fulfil these requirements, the distance of their vertex 
to the beam axis, R^, is calculated. The distribution of Ry is shown in Figure El for data from 
1994 and Monte Carlo: most of the photon conversions occur at radii between 40 and 90 mm, 
corresponding to the position of the two cylindrical layers of the SMD. Good agreement of the 



simulation of the photon conversion probabihty inside the TEC with the data is obtained after 
enlarging the conversion probability by a factor of about 1.6 for data taking periods after the 
installation of the SMD, to account for additional material not fully considered in the Monte 
Carlo simulation. The flat background stems mainly from 3-prong hadronic r decays. 

After rejection of identified photon conversions, 0.4% of the r decays still contain tracks 
from photon conversions. They are accounted for in the migration efficiencies determined from 
Monte Carlo. 

The effect of tracks scattered back from BGO clusters is investigated. As their momenta 
are low they are removed by the requirement on the transverse momentum of a track. 

Selection of e+e~ — > t+t~(7) events 

Events of the process e^e~ -^ r+r^(7) are characterised by two jets with low track and calori- 
metric cluster multiplicities, where a jet may consist of an isolated electron or muon. To ensure 
good track measurements only events in the barrel region of the detector are accepted by re- 
quiring I cos6thrust\ < 0.7. The event multiplicity is defined as the sum of the number of tracks 
in the TEC and the number of neutral calorimetric clusters, without an assigned charged track 
and with an energy larger than 0.5 GeV. This event multiplicity is required to be less than 10. 
Each event is divided into two hemispheres with respect to the plane orthogonal to the thrust 
axis. The main backgrounds arise from two-photon interactions, e"'"e~ events and Z decays into 
two muons. These processes are rejected using information from both hemispheres. 

• Two-photon interactions: each hemisphere must contain at least one calorimetric cluster. 
The sum of the energy deposited in the calorimeters and the muon momenta measured 
in the muon spectrometer must be larger than 13 GeV. 

• e"'"e~ events: the total energy deposited in the BGO must be less than 60 GeV. In addi- 
tion, the energy deposit in each hemisphere must be less than 44 GeV and the acoplanarity 
angle between the leading tracks of the two hemispheres must be larger than 0.003 rad. 

• Z decays into two muons: events with a track in the muon chambers with a momentum 
larger than 42 GeV are not accepted. Furthermore, events with a muon or a MIP in both 
hemispheres are rejected. 

A sample of 70016 e'^'e^ -^ r+r^(7) events is selected. The estimations of the efficiencies 
and background fractions are done separately for each year. The average selection efficiency 
for e+e~ -^ T^T~{'y) events inside the barrel, estimated from Monte Carlo, is 78.8 ± 0.2%. 

The background from hadronic and other leptonic Z decays and two-photon interactions is 
estimated from Monte Carlo. As an example, the distribution of the event multiplicity, is shown 
in FigureElfor data and a superposition of Monte Carlo from e^e" -^ T^T^{-f) and background. 
The background from hadronic Z decays dominates at large values of the multiplicity. After 
applying a correction factor to the fraction of the hadronic background of 1.05, very good 
agreement is obtained. The background from Bhabha events is determined using the upper part 
of the energy distribution measured in the BGO. This is shown in Figure|3]for an event sample 
with the Bhabha rejection cuts relaxed. The prediction from Bhabha Monte Carlo is scaled by 
a factor of 1.1 to agree with the data. The background after the final selection is estimated 
from Bhabha events using this scale factor. The same procedure is applied to estimate the 
//"•"/i" background. Figure shows the spectrum of muons measured in the muon spectrometer 



for data and Monte Carlo after applying a correction factor of 1.1. The background fractions 
from all sources are listed in Table H The cosmic background, estimated from the distribution 
of the distance of closest approach to the beam position, is found to be negligible. 

Determination of the topological branching fractions 

The maximum likelihood method is used to determine the topological branching fractions, with 
likelihood function L defined as: 

L = x{p{K,,, n:,^), (1) 

i=0 

where P is the Poisson distribution, A^^^,^ is the number of observed events and N^^^ is the 
number of expected events with i reconstructed tracks. The latter is: 

NLp = Nr E ^Oy^+E^S' (2) 

i=l,3,5 k 

where A',- is the number of r decays and B{j) is the branching fraction of j-prong r decays, 
j = 1, 3 or 5. The elements of the track detection efficiency matrix, e^^ , are determined from 
Monte Carlo. The non-diagonal elements represent migrations between the topologies. The 
number of non-tau background events, N^g, obtained from Monte Carlo, is normalised to the 
data luminosity. The index k runs over all background sources. 

The efficiency matrix of the track reconstruction is shown in TableEl indicating the numbers 
of reconstructed tracks for r decays to 1, 3 and 5 charged particles. Tracks arising from neutral 
kaon decays are not accounted for in the topology. Table El shows the number of observed 
r decays in the different topologies and the estimated background. In the fit the constraint 
B{1) + B{3) + B{5) = 1 is applied and the sum of N^.^^ is constrained to the number of observed 
r decays. The following results are obtained: 

B{1- prong) = 85.274 ± 0.105%, 
B(3 - prong) = 14.556 ± 0.105%, 
B{5- prong) = 0.170 ± 0.022%, 

where the uncertainty is statistical only. The x^/d.o.f. is 5.7/4. The correlation coefficients 
are given in Tabled 

In Figures El and 13 the number of observed r decays is shown as a function of the charged 
track multiplicity in linear and logarithmic scale, respectively. Also shown are the Monte Carlo 
predictions, using the fitted branching fractions, and the background. 

Systematic uncertainties 

The criteria to suppress the different background sources are varied within reasonable ranges 
and the changes in the topological branching fractions are taken as systematic uncertainties. 
The uncertainty on the cross section of e"'"e~ — *> hadrons |TB^ and the scale factor applied to the 
Monte Carlo normalisation are considered. The cross section uncertainties on e"'"e~ and ^i^ ii~ 
final states and two-photon interactions have negligible effects on the branching fractions. The 



background uncertainties from e"'"e~ and /x^/i^ final states are obtained from the statistical 
uncertainty of 1% on the scale factors applied to the Monte Carlo distributions shown in 
Figures ID and The systematic uncertainty due to track efficiency is obtained by varying 
this quantity within its statistical uncertainty of 0.25%. Furthermore, the track definition 
criteria are changed within reasonable ranges. The uncertainty from double track resolution 
is estimated by reweighting the A0 distribution of Figure Q forcing agreement between data 
and Monte Carlo in the low A0 region. The effects on the branching fractions are taken as 
systematic uncertainties. 

The systematic uncertainty due to photon conversions is obtained from the statistical un- 
certainty of 10% on the photon conversion probability correction factor and from variations of 
conversion identification criteria. 

The uncertainty from Monte Carlo statistics is included as statistical uncertainties on the 
efficiency matrix in Table 1. The resulting variations on the branching fractions are taken as 
systematic uncertainties. Effects from the energy scale uncertainties of the subdetectors are 
negligible. A summary of the systematic uncertainties is provided in Table El 

After combination of the systematic uncertainties the results for the branching fractions of 
the r lepton decays into one, three and five charged particle final states are: 

i3(r ^ (1 - prong)) = 85.274 ± 0.105 ± 0.073%, 
B{t ^ (3 - prong)) = 14.556 ± 0.105 ± 0.076%, 
B{t ^ {5 - prong)) = 0.170 ± 0.022 ± 0.026%, 

where the first uncertainty is statistical and the second is systematic. These new results are 
in agreement with a recent measurement with the full LEP statistics jS] and with the current 
world averages [T]. 
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Background source 


Fraction [%] 


Z -^ hadrons 


1.59 


Two-photon interactions 


0.16 


e+e~ — s> e"'"e~(7) 


0.16 


e+e" -^ /i+/i"(7) 


0.68 



Table 1: The background fractions from the different sources in the e"'"e 
sample. 



T ' T 



' (7) event 





^gen 


iVrec 


1 


3 


5 





7.58 ±0.02 


0.78 ±0.01 


0.84 ±0.15 


1 


70.18 ±0.05 


6.88 ±0.04 


3.58 ±0.30 


2 


0.33 ±0.01 


26.89 ±0.07 


9.55 ±0.50 


3 


0.16 ±0.01 


47.05 ±0.09 


20.50 ±0.73 


4 


<0.01 


0.21 ±0.01 


24.38 ± 0.79 


5 


<0.01 


0.05 ±0.01 


14.88 ±0.63 


6 


<0.01 


<0.01 


<0.01 



Table 2: The efficiency matrix of track reconstruction in percent. A^^en denotes the number of 
charged tracks of the r decay before detector simulation and N^-ec the number of tracks after 
the reconstruction. 



Nrec 





1 


2 


3 


4 


5 6 


Data 


11935 


107283 


8166 


12378 


216 


53 1 


Background [%] 


3.7 


1.9 


9.7 


2.5 


33.8 


7.6 



Table 3: The number of r decays observed in the different topologies and background estimated 
by Monte Carlo. 



Source 


B{1 — prong) 


B{3 — prong) 


B{5 — prong) 


B{\ — prong) 


1.0 


-0.978 


-0.082 


B{3 — prong) 




1.0 


-0.127 


B{5 — prong) 






1.0 



Table 4: The correlation coefficients obtained from a ffi of the topological branching fractions. 



10 



Source 


B{1 — prong) 


B{3 — prong) 


B{5 — prong) 


Z -^ hadrons 


0.048 


0.052 


0.024 


e^e~ -^ e'^e~(7) 


0.010 


0.010 


0.001 


e+e" ^ fi+fi-{-f) 


0.010 


0.010 


0.001 


Two-photon interactions 


0.011 


0.011 


0.001 


Track definition 


0.035 


0.035 


0.003 


Double track resolution 


0.012 


0.012 


0.001 


Photon conversions 


0.017 


0.017 


0.004 


Monte Carlo statistics 


0.032 


0.032 


0.007 


Total 


0.073 


0.076 


0.026 



Table 5: Systematic uncertainties in % on the branching fractions resulting from the listed 
sources and their combined values. 
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Figure 1: Distribution of the azimuthal angle between two adjacent tracks A0, where the 
tracks must have a transverse momentum larger than 10 GeV. 
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Figure 2: Distribution of the radial distance from the beam axis, i?„, of vertices reconstructed 
using photon conversion tracks. The flat hatched distribution stems from track pairs of hadronic 
r decays. 
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Figure 3: The distribution of the event muhiphcity. The Monte Carlo prediction for 
e"'"e~ -^ r^T^{'-f) and the background from other leptonic and hadronic Z decays after ap- 
phcation of the scale factors is also given. 
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Figure 4: The distribution of the BGO energy for selected e"'"e~ -^ r"'"r~(7) events with relaxed 
cuts against e"'"e~ -^ e^e~ (7) background. Also shown is the Monte Carlo expectation for 
e+e" — > T^T~{^) and e"'"e~ — > e"'"e~(7) events after rescaling. 
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Figure 5: The momentum distribution of tracks measured in the muon chambers in the 
e"'"e~ -^ r^T^{'-f) event sample with relaxed cuts against dimuon background. Also shown 
is the Monte Carlo expectation for e^e~ -^ r^T^{-f) and the background from dimuon final 
states after rescaling. 
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Figure 6: The charged multiphcity distribution from r decays. Also shown is the expectation 
from Monte Carlo for e^e~ — > T^T~{-y) and the background. 
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Figure 7: The charged multiphcity distribution from r decays. Also shown is the expectation 
from Monte Carlo for e"''e- -^ r+r-(7) and the background from hadronic Z decays and other 
sources. 
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